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^ Abstract 
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We present a scheme to couple trapped 87 Rb atoms to a superconducting flux qubit through a 

a 

magnetic dipole transition. We plan to trap atoms on the evanescent wave outside an ultrathin 
fiber to bring the atoms to less than 10 urn above the surface of the superconductor. This hybrid 
setup lends itself to probing sources of decoherence in superconducting qubits. Our current plan 
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CT) has the intermediate goal of coupling the atoms to a superconducting LC resonator. 
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INTRODUCTION 



A hybrid quantum processor exploits the best aspects of its constituent qubits: atomic 
systems exhibit excellent isolation from the environment, and condensed matter systems 
benefit from technological developments in microfabrication, leading to a scalable qubit 
architecture pp. No known individual quantum system satisfies all the DiVincenzo criteria for 
a quantum computer [2]; meeting these requirements is the driving force behind developing 
hybrid systems. 

Here we propose a scheme to couple the ground state hyperfine transition of neutral 87 Rb 
atoms to a superconducting (SC) flux qubit. A flux qubit consists of three or more Josephson 
junctions in a micrometer-sized superconducting loop[3H3]. The supercurrent forms a super- 
position of clockwise and anticlockwise flow. Typical designs involve one junction smaller 
than the others so that its Josephson inductance dominates the loop inductance. We can 
use the microwave magnetic field from this SQUID loop to drive the atomic transition. The 
flux qubit possesses a long coherence time, on the order of 10 /is [6] at the flux degeneracy 
point in a four-junction SQUID [TJ, and a tunability range of over a GHzjB]. 

The atomic qubit consists of two magnetic sublevels on one of the hyperfine ground state 
manifolds of 87 Rb, between the \5Si/ 2 ; F = 1) and 155*1/2; F = 2) states, which are separated 
by about 6.83 GHz. 

Similar proposals focus on creating a memory of ions [JJ H3], neutral atoms [TUl ITT] , or 
molecules [T2"HH] to couple to a SC stripline resonator [T5]. Nitrogen- vacancy (NV) centers 
have taken a prominent role in these designs because of the ease with which these natural 
ion traps integrate with SC systems [16]. Recent results demonstrate strong coupling of an 
ensemble of NV centers to a stripline resonator [17] . 

THE INTERACTION AND ITS REALIZATION 

The flux qubit induces magnetic dipole transitions in the ground state hyperfine manifold 
of 8 Rb. The coupling occurs via the interaction Hamiltonian, ^ — fii ■ B, where Hi is the 
magnetic dipole moment of the zth individual atom, and B is the amplitude of the microwave 
magnetic field from the SC qubit at the same atomic frequency. To calculate the strength 
of the coupling we must first determine the strength of the magnetic field. 
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The average magnetic field associated with a single microwave photon is: 

where u is the frequency of the photon, /i the permeability of free space, and V e g the effective 
mode volume. We take the effective mode volume for a square SQUID with 10 /zm x 10 /im 
sides and width of 5 fim to be about 10 x 10~ 16 m 3 , assuming the field is confined within 
5 micrometers above and below the strip. For a photon at 6.8 GHz the magnetic field is 
about lxlO -8 T, which for the typical moment of 1.4 x 10 10 Hz/T gives a coupling strength 
of roughly 100 Hz for one atom. 

The field associated with a single quantum flux, $ = |^ = 2 * 10~ 15 Tm 2 , in a similar 
device corresponds to 2 * 10 _5 T, a much higher value than the single photon. The coupling 
of one magnetic flux to the atomic magnetic moment would be lowered by some geometric 
factor determined by the specific shape of the magnetic field obtaining a single flux coupling 
larger than 100 Hz. Further understanding of the field distribution and the inductance of 
the SQUID will help narrow the range of this number. 

We are working to first couple atoms to a high-Q lumped-element superconducting LC 
resonator, tuned to the 87 Rb hyperfine splitting. Figure 1 presents our current lumped- 
element LC resonator. The niobium resonator consists of a meandering inductor and an 
interdigitated capacitor coupled to a transmission line. At a working temperature of 12 mK 
and on resonance at 6.863 GHz, the transmission through the microwave line decreases by 
1.5 dB, and the loaded quality factor is 40,000. We simulate the electromagnetic fields of 
the SC resonator with the software package, High Frequency Structure Simulator (HFSS), 
to determine the position of maximal magnetic field and uniformity. The coupling to the 
resonator depends on the specific geometry but will not reach the single photon value that 
we estimate for the SQUID. It will be lower by at least an order of magnitude. Figure 2 
shows that the magnetic field outside of the SC as the field lines encircle the edges of the 
meandering inductor. Using HFSS, we integrate the electromagnetic field over the entire 
simulation region to find the total energy with a total photon number of n p h— 0.016 photons. 
Then scaling the magnetic field located 5 /im above the inductor by 1/ y/n p h we obtain the 
magnetic field produced by one photon as 2.47 • 10~ 9 T. We expect the coupling per atom 
to be 34.6 Hz. 

The resonator operates inside a Triton 200 Cryofree Dilution Refrigerator from Oxford 
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Instruments at 12 mK. The 200 fiW cooling power of the final stage and degradation of 
the quality factor of SC resonators in a strong magnetic field present an unusual set of 
requirements for trapping atoms. These design constraints make the magneto-optical trap 
(MOT) and dipole trap difficult to implement [HJ EE]. We will realize an evanescent wave- 
based dipole trap outside an optical fiber whose diameter is smaller than the wavelength of 
input light. This trap is now called a tapered optical fiber (TOF) trap [20TT22] (see Fig. 3). 

THE SUPERCONDUCTING RESONATOR 
Tuning Scheme 

We have developed a frequency tuning system for a "lumped-element" thin-film supercon- 
ducting Nb microwave resonator on sapphire for coupling to cold 87 Rb atoms. The resonator 
must be tuned in order to match the resonance to the energy splitting of the ground state 
of 87 Rb. We employ an Al pin as a frequency tuner by placing it above the inductor and 
using a piezo-electrically driven stage to change the inductance of the resonator. 

Tuning Performance 

Figure 4 displays the data for tuning the resonator. The present resonator has a large 
meandering inductor to ease alignment issues. 

The data shows a decrease in the frequency as the tuning arm approaches the resonator 
(see Fig. 4). The first design allowed the pin to float as the piezo stage moves. The 
ungrounded pin showed an unexpected shift of 66 MHz over a 290 fim range of motion, 
while lowering the loaded Q from 87,000 to 18,000. Connecting a copper foil from the pin 
to the sample box rectifies this grounding issue. 

EXPERIMENTAL APPARATUS 

Figure 5 shows the current apparatus to test the LC resonator and effects of optical fibers 
on the resonator. It consists of a copper sample box thermally anchored to the bottom 12 
mK stage of the refrigerator. 
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It contains a piezo-electric driven stage to adjust the position of a tuning wire above the 
LC resonator. We are currently designing the inital cooling scheme necessary to load atoms 
into the TOF trap. The ultrathin fiber would run through the cooled atoms allowing for 
their transfer to the TOF trap. The TOF trap introduces the atoms into the sample box 
and is located 10 /xm above the resonator providing for a coupling of about 10 Hz per atom. 

THE TAPERED OPTICAL FIBERS 

Intensity Profile and Trapping Potential 

Figure 6 displays the intensity profile of the fundamental mode of a linearly y-polarized, 
980 nm input laser about a 500 nm diameter fiber. Notice the discontinuity of the field when 
it leaks outside the glass fiber. The intensity profiles inside and outside the fiber are [231 El] 

\E m \ 2 = g m [J 2 {hr) + uJl{hr) + fj 2 (hr) 

+ (uJ 2 (hr) - f p J (hr)J 2 (hr)) cos [2 (0 - O )]] (2) 
\E out \ 2 = g out [K 2 (qr) + wK 2 (qr) + jK\{qr) 

- (wK 2 (qr) - f p K (qr)K 2 (qr)) cos [2 (0 - O )]] , (3) 

where J n and K n are Bessel functions of the first and second kinds of order n respectively, 
r is the radius from the center of the fiber, and g is a normalization constant. The terms 
K q, fi fp and w, given explicitly in [23], are functions dependent on the fiber radius and the 
propagation constant, (3, of the input field. 

The intensity outside the fiber falls off exponentially with a characteristic decay length 
proportional to the wavelength of the laser light (see Fig. 6). We can realize a trapping 
potential by sending two color light, red- and blue-detuned from the D 2 transition of Rb, 
through the fiber [20] • The corresponding attractive and repulsive potentials combine to 
create a trap depth on the order of mK. 

Using linearly-polarized trapping light and coupling the red-detuned laser from both ends 
of the fiber, we will form ID optical lattices along either side of the fiber waist. The lattice 
sites admit either zero or one atoms due to collisional blockade effects, yielding ~10,000 
87 Rb atoms within a few hundred nm from the waist of a 500 nm TOF [25]. Figure 7 shows 
a calculation of the trapping potential for a fiber with a 500 nm diameter. We apply lasers 
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with wavelengths of 730 and 980 nm, and powers of 30 and 13 mW, respectively. These 
parameters produce a 

trap depth of 7.46 mK with the atoms trapped 137 nm from the fiber surface. The ratio of 
red- to blue-detuned light controls the trap depth and position of the minimum. Increasing 
the red-detuned power, deepens the trap, and the atoms approach the fiber surface. 

The proximity of the atoms to the fiber surface means we must include the van der Waals 
interaction in the trapping potential. We do not, however, use the exact van der Waals 
interaction between an atom and a nanowire (26] but instead treat the system as if the 
atoms were located next to an infinite dielectric, yielding a C3 coefficient of 8.46xlO _49 J- 
m 3 [27] 128]. Given the distances of atoms from the fiber surface, it appears that we may be 
in a cross-over region between the van der Waals and Casimir-Polder regimes. To reconcile 
this issue we compare the asymptotic solution of the Casimir and Polder interaction [2H1 EH] , 
which accounts for retardation effects, to the dipole-dipole van der Waals term. 

Figure 8 displays the van der Waals and Casimir-Polder potentials. For this calculation 
we assumed a value of 0.0794/i Hz-cm 2 /V 2 for the ground state polarizability of Rb, and of 
2.04 for the dielectric constant of silica, where h is Plank's constant. Using Casimir-Polder, 
rather than van der Waals, we find a trap depth of 7.43 mK with intensity minima located 
at 138 nm. This results in a shift in trap depth of 0.03 mK and atomic trapping sites of 1 
nm, so we can neglect the Casimir-Polder interaction. 

Advances on the Fiber and its Realization 

The fiber pulling device necessary to construct a TOF requires a hydrogen-oxygen flame 
to create a local heat zone, while stepper motors pull the ends of a normal fiber until 
achieving the desired fiber waist radius. The TOF consists of three regions: the unmodified 
fiber, the taper or transition, and the fiber waist [21] ■ How one modifies the length of the 
hot zone during the pulling process governs the shape of the TOF [32J. The geometry of 
the transition region affects how light couples between the unmodified region and the fiber 
waist. Poor coupling leads to scattered light, which heats the SC. The design of the taper 
must satisfy an adiabatic criterion that best transfers light from the propagating mode in 
the unmodified region, LP 01 , to the HE n mode in the fiber waist [33]. Reference [31] gives 
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the adiabatic condition as 

n < P(h-h\ ( 4) 

where Q is the angle of the taper, p the radius of the fiber along the taper, and /3i j2 the 
propagation constants of the fundamental and first excited modes, respectively. When the 
light propagates in the unmodified region the core-cladding interface guides the mode, and 
the core contains a majority of the light. As the fiber tapers, the core thins, and eventually 
the mode spreads into the cladding. This radius change leads to a breaking of the single mode 
condition, and higher modes may be excited. The radius of the transition region continues to 
decrease until the cladding- air interface guides the light. When the cladding radius reaches 
a size that admits only a single mode, this acts as a mode filter, and all the higher modes 
scatter or reflect. To avoid leaking light and heating the superconducting surfaces, the 
tapered region of the fiber must have a shape that satisfies the adiabatic condition of Eq. 4. 

OPEN QUESTIONS 

We are still investigating options to introduce the atoms to the TOF trap in the dilution 
refrigerator. Currently we plan to anchor a hollow copper rod to the 4K stage in the dilution 
refrigerator placed beside the sample box. Inside the rod we will explore one of three options: 
a pyramid MOT [35], a diffraction grating-based pyramid MOT [36], or a trap relying on 
optical pumping [37]. The first two MOTs provide a compact design and require only a 
single laser. The final option offers a trap without magnetic fields. An electron beam driven 
atomic source, which requires only 10 mW of power [3E], can introduce the atoms to the 
MOT. The fiber waist runs through the center of the trap allowing for the transfer of atoms 
from the MOT to the TOF trap where an optical conveyor belt can move the atoms above 
the flux qubit or the LC resonator [39|. 

CONCLUSION 

We propose a scheme to couple 87 Rb atoms, trapped along a TOF, to a flux qubit via the 
magnetic dipole moment. As an intermediate step the atoms will couple to a SC lumped 
element resonator. This project has applications in quantum information by working towards 
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a hybrid quantum computer. Trapping atoms near a superconducting circuit also presents 
an opportunity to characterize sources of charge noise in the decoherence of SC qubits. 
Specifically, we can probe the source of charge noise by exploiting the atoms, or use Rydberg 
atoms, with their large dipole moments, as electric field sensors. This system can help to 
uncover the microscopic agents responsible for decoherence in SCs and use this to increase 
the coherence time. 

Work supported by National Science Foundation of the USA through the Physics Frontier 
Center at the Joint Quantum Institute. 
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FIG. 1. (Color online) Current lumped-element LC resonator. 1: microwave transmission line. 2 
interdigitated capacitor. 3: meandering inductor. 
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FIG. 2. (Color online) HFSS simulation of the H field produced by 0.016 photons in the lumped- 
element LC resonator. 
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FIG. 3. (Color online) A two color TOF trap hold atoms 10 /jm above a flux qubit (small green 
square on the superconducting chip). 
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FIG. 4. (Color online) The shift in the resonance frequency of the lumped-element resonator as a 
function of how far the aluminum pin moves. The blue, red, green, and black lines represent the 
aluminum tuning arm moving 0, 0.1, 0.2, and 290 //m, respectively. These shifts yield resonances of 
6.8637, 6.8618, and 6.8571 GHz with loaded Q's of 87,000, 24,000, 19,600, and 18,700, respectively. 
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FIG. 6. (Color online) Calculation of the intensity profile for a 500 nm fiber with an input of 980 
nm linearly polarized light. 
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FIG. 7. (Color online) The trapping potential for red- (long-short dashed line) and blue-detuned 
(dashed line) light of 980 and 730 nm with input powers of 30 mW and 13 mW, respectively. The 
short dashed line represents the van der Waals interaction and the solid line is the total combined 
potential. These parameters yield a trap depth of 7.46 mK, 137 nm from the fiber surface. 
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FIG. 8. The solid line represents the van der Waals potential while the dashed line denotes the 
C asimir- Polder . 
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